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Abstract 



Using the axial- vector coupling and the electromagnetic form factors of the D and D* mesons in 
4^ 2+1 flavor Lattice QCD, we compute the D*D7r, DDp and D*D*p coupling constants, which play 

^ an important role in describing the charm hadron interactions in terms of meson-exchange models. 

We also extract the charge radii of D and D* mesons and determine the contributions of the light 
and charm quarks separately. 
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I. INTRODUCTION 



An approach that utihzes the effective meson Lagrangian can successfully describe the 
suppression of J /ip production, which is considered to be a signal for the formation of quark- 
gluon plasma in relativistic heavy-ion collisions (RHIC) [IHl] . This model has also been used 
to study J /ip absorption by nucleons [5], charm production from proton-proton collisions |6] 
and charm photoproduction off the nucleons |7]. The effective Lagrangian employed in this 
method includes interaction vertices among tt, p, J/'?/', D and D* . The coupling constants 
of these mesons then play a crucial role in giving an accurate description of charm-hadron 
production and suppression in collisions performed at RHIC. 

In constructing the effective Lagrangian, meson exchange models use SU(4) symmetry 
together with Vector Meson Dominance model (VMD) to determine various coupling con- 
stants. Since the coupling of D and D* to the pion occurs only through anomalous vertices, 
the dominant contribution to processes involving diagonal transitions of D and D* mesons 
comes from p exchange. Our primary aim here is to calculate the DDp and D*D*p coupling 
constants from first principles using 2+1 flavor Lattice QCD. As a byproduct, we obtain the 
electromagnetic form factor and the charge radii of D and D* mesons. In order to benchmark 
our simulations and compare with the available literature, we started our calculations with 
the well-known D* Dtt coupling constant, gD*D-K- Therefore we give our results for gD*Dir 
also for illustrational purposes. 

The calculation in this work is reminiscent of the pion (or kaon) electromagnetic form 
factor, which is considered to be a good observable to test QCD in a broad range of energy 
regime. There are also similarities between the p (or K*) [8] and the D* electromagnetic 
form factors considered here. It has been found that the experimental data of the pion 
electromagnetic form factor at small momentum transfer are described quite successfully by 
the VMD ansatz P 

where rUp is the p-meson mass, (^ttttp and Qp are the vr-p and p-photon coupling constants, 
respectively. While the p-meson is expected to dominate the electromagnetic current around 
the meson pole in the timelike region, VMD can describe the data quite accurately up to 
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~1 GeV^ in the spacelike region [10]. The monopole form, 

i^mon(Q') = i^g2/A2 (2) 

inspired by the VMD hypothesis (by assuming universahty (^TrTrp = fl'p), has also been used 
as a useful tool to fit the data to and predict the charge radius of the pion p!THT3] . It has 
been inferred from a compilation of experimental and theoretical results that the deviation 
from VMD starts above ^1.5 GeV^ [15], which can be identified as the transition energy 
scale for pion from low-energy behavior to perturbative QCD. Since the D meson is much 
heavier than the pion, the transition is expected to occur at higher momentum transfers. 
Motivated by the success of VMD in describing the electromagnetic form factor of the pion, 
we use the same method to calculate the DDp and D*D*p coupling constants from lattice 
QCD data. 

0*071 coupling constant has been determined experimentally by CLEO Collaboration as 
Qd'Dtt = 17.9 ± 0.3 ± 1.9 and studied in the literature extensively. Therefore this observable 
can serve as a useful benchmark tool in this sector. While the results for gD*Dn from early 
QCD sum rules studies are well below the experimental value [151 IIS]) those from lattice- 
QCD works are in good agreement with the experiment [T71[TS]. Our calculations for go'Dn 
here improve upon previous studies in several aspects, such as the lattice size and the number 
of sea-quark fiavors. 

Our work is organized as follows: In Section [IT] we present the theoretical formalism of D 
and D* form factors together with the lattice techniques we have employed to extract them. 



In Section III we give and discuss our numerical results. Section |IV| contains a summary of 



our findings. 

II. THE FORMULATION AND THE LATTICE SIMULATIONS 

We compute the meson matrix elements of the vector current = |c7^c+|m7^u — l(i7^(i, 
which can be written in the form 

{D{p')\V,{q)\D{p)) = ip + p'), [e^F'^iQ') + e,F%Q')] (3) 
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for the D meson. As for the spin-l D* meson, we have 
{D*ip',s')\V,iq)\D*iP,s)) = e':ip',s') \G^iQ')g^'^{p^ +pn 



(4) 

where e and e' are the polarization vectors of the initial and final vector mesons, respectively. 
The form factors Gi 2,3 can be arranged in terms of Sachs electric, magnetic and quadrupole 
form factors as follows 1191: 



Fcm = G,m + l^FQm 

Fm{Q') = (5) 
= G,{Q') - G,{Q') + (1 + v)G,{Q') 

where rj = (5^/4m|,,. 

Here we consider and D*"*" mesons therefore we take Cq = 1/3 (anti-(i-quark) and 
Cc = 2/3 (c-quark). We use the notation = —Q^ = {p — p'Y ^-^id, F'^ and F"^ are the vector 
form factor of where the external field couples to the c- and the d-quark in the D-meson 
respectively. In the limit of vanishing four-momentum transfer, we have F'^(O) = F''(0) = 1. 
Note that -Fem(Q^) = [^cF'^iQ'^) + ^qF'^iQ'^)] gives the electromagnetic form factor of D and 
we have -Fem(O) = 1 due to charge conservation. Similar constraints hold also for the D* 
electric form factor, F^j^ = Fc(Q^). 

The D*D7i coupling constant, go'Dn, can be accessed via the transition matrix element 
{D{p')\A'^{q)\D*{p, s)), where the axial-vector current for the light quark is given by = 
M757^(i. This matrix element can be parameterized with three form factors, Fo(g^), Fi{q^) 
and F2{q^): 

{D{p')\A>^{q)\D\p,s)) = 2mvFo{q'f-^q^ 

+ (mn + mn*)F^iq^)[e''^ - '^q^] (6) 

+ F,{q^) ^^-^ [p^+pi^ - 

rriD + mD* q^ 

PC AC relation and the VMD imply that the divergence of the axial- vector current q^^A^ is 
dominated by a soft pion: 

{D{p')\q^A^{q)\D*{p,s)) = g^.^^^^M±^f^^l + . . . . (7) 



We refer the reader to Ref. [T7] for the details of our calculations for Qwd-k] here we just 
summarize the main steps. Rewriting the matrix element in terms of transferred and final 
momenta by defining p^^ = {p' + g)^, we can identify go'Dn in terms of Fi(0) and -^2(0) as 



Defining 



go^Drr = Y i(^D + m^.)^i(O) + {rriD* - mfl)F2(0)] . (8) 

JIT 



Jtt Jtt 



and rearranging Eq. ([8| we write Qd'O-k as 



9D*D. = GM[^ + ^^,]- (10) 



(11) 



Gi(o); ■ 

To extract the coupling constants, we compute the mesonic two-point 
(C(t;p)) = 5^e-P-''(vac|T[x(a;)x(0)]|vac), 

X 

{C^uit-v)) = 5^e-^P-^(vac|r[x^(x)x.(0)]|vac), 

X 

and three-point correlation functions, 

((7°(t2,ti;p',p)) = E e-*P-^^e^'^-^Mvac|T[x(x2)\/"(xi)x(0)]|vac), (12) 

X2,Xl 

(C'^^(t2,ti;p',p)) = E e-^P"'^e*''-^Mvac|T[x^(x2)V"(xi)x.(0)]|vac), (13) 

X2,Xi 

(C'^.(t2,ti;p',p)) = E e-^^-''^e^''-''Mvac|T[x(x2)A^(xi)x.(0)]|vac). (14) 

X2,Xl 

The meson interpolating fields are given as 

X{x) = [c/(x)75c(x)], x^iix) = [c/(x)7^c(x)]. (15) 

In our setup, the three momentum of the outgoing meson is automatically projected to zero 
momentum due to wall method, which is explained below; i.e. p ' = 0. 

In terms of the quark propagators S{x,x'), the three-point correlator for the D meson 
(we take the vector-field coupling as an example) can also be written as 

{C^{h,h;p',p)) = -^J2 e-'^'^^e''^--- (16) 

X2,Xl 

x(Tr[75 ^d(0,xi) 7" Sd{xi,X2) 75 S,{x2M)- (17) 
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A similar expression holds also for the D* meson. While point-to-all propagators Sd{0,xi) 
and Sc{x2,0) can be easily obtained, the computation of all-to-all propagator Sd{xi,X2) is 
a formidable task. One common method is to use a sequential source composed of Sd{0,Xi) 
and Sc{x2, 0) for the Dirac matrix and invert it in order to compute Sd{xi, X2) [20\. However, 
this method requires to fix sink operators before matrix inversions. 

An approach that does not require to fix sink operators in advance is the wall method, 
where a summation over the spatial sites at the sink time point, X2, is made before the 
inversion. This corresponds to having a wall source or sink: 

(Q^(t2,ti;0,p)) = -2 Yl ^"""^ 

x(Tr[75 Sd{0,x,) r Sd{x,,x'2) 75 S,{x2,0)]) 

where the propagator (instead of the hadron state) is projected on to definite momentum 
{S and W are smearing labels for shell and wall). Since the wall sink/source is a gauge- 
dependent object, one has to fix the gauge. Here we fix the gauge to Coulomb which produces 
a somewhat better coupling to the hadron ground state as compared to the Landau gauge. 
The wall method has the advantage that one can first compute the shell and wall propagators 
and then contract the propagators in order to obtain the three-point correlator, avoiding 
any sequential inversions. Use of the wall method allows us to compute the D and D* axial- 
vector and electromagnetic-transition channels simultaneously, which would require separate 
treatments with the traditional sequential-source method. 
We compute the matrix element in Eq. ([s]) using the ratio 



^"(i2,ti;p ,p;/i) = 

(C'g^y(t2,tl;p^p)) 
{Csw{t2;p')) 



{Cssit2 - ti; p)) {Cssjti; p')) {Css{t2; pQ) 
{Css{t2 - ti; pO) {Cssih; p)) (^55(^2; p)) 



1 1/2 (19) 



ti is the time when the vector field interacts with a quark and ^2 is the time when the final 



meson state is annihilated. The ratio in Eq. ( 19 ) reduces to the desired form when t2 — ti 
and ti ^ a, viz. 

where m£, and are the mass and the energy of the initial baryon. We apply a procedure 



of seeking plateau regions as a function of ti in the ratio (20 ) and calculating the vector form 
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factors F%Q'^) and ^'(Q^). We extract the D -meson mass from the two-point correlator 
with shell source and point sink, and use the dispersion relation to calculate the energy at 
each momentum transfer. 

The matrix element in Eq. ^ is computed using the ratio 



{C'swit2;p')) 



This ratio gives 



'33 



(gg5fe-ti;p)) (gg-(ti;pO)(Cg-(t,;pO) 

(C^5(^2-ti;pO) (C^5(^i;p))(C^5(^2;p)) 



1/2 



(21) 



In order to single out the electric form factor we compute 



(22) 
(23) 

(24) 



Finally, we compute the form factors -Fi(O) and ^2(0) needed for gD*D-w via the ratios [T7] 



Ri{t) 



Rz{t) 



R2{t) 



Ri{t) 



(25) 



(fjCwwih — ti) ' C'||(ti, qjCwwih — ti) 

The masses and the normalization factors Z^* and are obtained from exponential fits to 
the zero-momentum two-point correlators. 



C^^{U-p)^Zd* 



-En*ti 



(26) 



Fi{0) can be computed easily, however we should compute -^2(0) by extrapolation since the 
term including F2{q^) vanishes at zero momentum transfer: 



Fi{q = 0, U] 



Ri{t 



rriD' + rriD 



FAq 



2tx 



Riiti 



mi:,* + rriD 



[mp* + mp) 
2m|,g*^ 



[q - Ed,{Ed, - mD>) ^ t^^t^ + i 



Ed* 



Ri{t\ 



Ed* Ri{ti) 
(27) 
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We assume the value of F2{0) to be close to its value at the smallest finite momentum transfer 
since we expect the F2/F1 ratio to be insensitive to the changes of transferred momentum 
around the pion pole. 

We make our simulations on a 32'^ x 64 lattice with 2+1 flavors of dynamical quarks and 
the gauge conflgurations we use have been generated by the PACS-CS collaboration [21] with 
the nonperturbatively 0(a)-improved Wilson quark action and the Iwasaki gauge action. We 
use the gauge conflgurations at /3 = 1.90 with the clover coefficient csw = 1-715, which give a 
lattice spacing of a = 0.0907(13) fm (a^^ = 2.176(31) GeV). The simulations are carried out 
with four different hopping parameters for the sea and the u,d valence quarks, Hgea, ~ 
0.13700, 0.13727, 0.13754 and 0.13770, which correspond to pion masses of ~ 700, 570, 410, 
and 300 MeV. The hopping parameter for the s sea quark is flxed to k^^^ = 0.1364. 

It is well known that the Clover action has discretization errors of 0{mqa). Precision 
calculations such as the spectral properties and the hyperflne splittings may require removal 
or at least suppression of these lattice artefacts by considering improved actions such as 
Fermilab [22] • On the other hand, calculations which are insensitive to a change of charm- 
quark mass are less demanding in this respect [23]. Considering also the precision levels 
we aim for the coupling constants and the flne spacing of our lattice, we choose to employ 
Clover action for the charm quark. We have checked the variation in our results by changing 
the charm-quark mass mildly and we conflrmed that the coupling constants are insensitive 
to such a change (see the discussion below). Note that the Clover action we are employing 
here is a special case of the Fermilab heavy quark action with csw = ce = cb [21] ■ We 
determine the hopping parameter of the charm quark {kc = 0.1224) so as to reproduce the 
mass of J/ip. 

We employ smeared source and wall sink, which are separated by 12 lattice 
units in the temporal direction. Source operators are smeared in a gauge-invariant 
manner with the root mean square radius of ~ 0.5 fm. All the statistical er- 
rors are estimated via the jackknife analysis. We make our measurements on 
45, 50, 50 and 70 conflgurations, respectively, for each quark mass. For the 
D-meson vector coupling, we make nine momentum insertions: {\Px\, \Py\, \Pz\) = 
(0, 0, 0), (1, 0, 0), (1, 1, 0), (1, 1, 1), (2, 0, 0), (2, 1, 0), (2, 1, 1), (2, 2, 0), (2, 2, 1). Given Eq. 
for the D* meson, we are limited with the kinematics where the momentum should be 
inserted in one direction only. For the D* meson, we make our measurements with four 
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momentum insertions: (bxl, byl, bzl) = (0, 0, 0), (1, 0, 0), (2, 0, 0), (3, 0, 0). We then rotate 
momentum in other directions and using the isotropy of space we average over equivalent 
momenta for both D and D* in order to increase the statistics. 
We consider point-spht lattice vector current 

= l/2[q{x + /i)f/t(l + 7^)g(x) - q{x)U^{l - ^,)q{x + /i)], (28) 

which is conserved by Wilson fermions. Therefore it does not require any renormalization on 
the lattice. The local axial-vector current, on the other hand, needs to be renormalized on 
the lattice, where the renormalization factors are computed in a perturbative manner [25] . 



III. RESULTS AND DISCUSSION 



We begin our discussion with the D*Dti coupling constant. We find that the dominant 
contribution to the coupling constant comes from the first term in Eq. ([s]). The second 
term including the ratio G2/G1 contributes to the coupling around 10%. Our results are 
listed in Table [Tj We give the dominant and minor contributions to the coupling constant 
individually. 





Gi(g2 = 0) GsM QD'D. 


0.13700 


14.15(1.58) 0.09(2) 15.45(1.78) 


0.13727 


13.63(1.57) 0.12(4) 15.24(1.81) 


0.13754 


12.76(1.43) 0.15(7) 15.54(2.08) 


0.13770 


15.46(2.17) 0.07(6) 16.44(2.41) 


Lin. Fit 


16.23(1.71) 


Quad. Fit 


17.09(3.23) 



TABLE I: Dominant and minor contributions to the coupling gD*Dn at each sea-quark mass we 
consider. 



As for the vector couplings, we construct the ratios in Eqs. (19) and (21) in order to 
extract the form factors. To demonstrate the signal-to-noise ratio and show the chosen 
plateau regions, we give in Fig. [ijthe ratio in Eq. (19) (Eq.(21)) as functions of the current 



insertion time, ti, for k = 0.13700 and for first seven (four) momentum insertions in the case 



9 




FIG. 1: The ratio in Eq. (19) (Eq.(21)) as functions of the current insertion time, ti, for k 
0.13700 and first seven (four) momentum insertions in the case of D (D*). 



To extract the coupling constant to the p-meson we use the VMD approach. There are 
different versions of VMD (see Ref. for a review). The two versions, which were first 
discussed by Sakurai [9J, differ by the mechanism the photon interacts with the hadron. In 
the more popular version, the photon is not allowed to directly couple to the hadron but 
only through a p-meson. This yields the following expression for the vector form factor: 

,2 



FviQ' 



QODp 



(29) 



where Qp is a constant which determines the coupling of the vector meson to the photon. In 
this version in order to satisfy the constraint -F(O) = 1 one has to to assume goDp = Qp- In 
a second version of VMD, the photon can couple to both the hadron and the p meson: 



FviQ' 



1 - 



9DDp 

ml + Q'^ Qp 



(30) 



In this version F(0) = 1 is automatically satisfied and one does not need to assume goDp 



Qp. We shall use the form in Eq. (30) to fit our data to and extract the coupling constants 



goDp and gD*D*p- We note that at the meson pole — )■ — m^, both the VMD forms are in 
agreement giving m^/(Q^ + m^)^^^^, so that the difference comes (if goDp 7^ Qp) at — )■ 



extrapolation. 
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The coupling constant gp can be obtained from partial decay width of the p-meson to 

e~^e~ , 

T{V ^ e+e-) = ^a^^ (31) 

with the fine structure constant a = 1/137. Using the experimental information [27j we find 
Qp = 4.96. We shall neglect the contributions from the p — u mixing as we have exact isospin 
symmetry on our lattice and such contributions are expected to play a role in the space-like 
region around the p-meson pole. 

We can obtain the electromagnetic charge radius of the D and D* mesons from the slope 
of the form factor at = 0, 



(32) 



For the monopole form in Eq. (|2| we have 



(-^) = ^- (33) 
Inserting this expression back into Eq. ([2]) for and rearranging we obtain 

{'•^) = 4 - l) . (34) 



<3" V / 

We extract the charge radii of the D and D* mesons using the above expression at the 
lowest finite momentum. This is a similar approach to that used in Ref. |8] to calculate 
light meson electromagnetic form factors. We calculate the light-quark and charm-quark 
contributions to the charge radii separately. This information helps us to study the charge 
radii of individual quarks that are bound in the meson. 

In Fig. [2] we show our lattice data for the form factors F^*^'^{Q'^)^ F^*'^^{Q'^) and the 
electromagnetic form factor -Fem(Q^) functions of the four-momentum transfer for 
i^u,d = 0.1370 with their jackknife errors. While we show data up to ~ 1.50 GeV^, we 
make our fits to first seven momentum-transfer values only (up to 1 GeV^). This is a 
region where VMD is expected to be valid. Nevertheless, the data at higher momentum 
transfers are also well described by the VMD form (for F'^*^'^{Q'^)) and the monopole form 
(for the electromagnetic form factor) as can be seen in Fig. |2} 

Our numerical results are provided in Table [11} We give the values of the coupling con- 
stants goDp and gD*D*p at four different light quark masses and as determined by a fit to the 
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FIG. 2: The form factors F^*^'^{Q'^), and the electromagnetic form factor Fem(<3^) as a 

functfon of the four-momentum transfer (in lattice units) for k = 0.1370. The shaded regions 
show jackknife error bars. 



VMD form in Eq. (30). We also give the values of the charge radii of D and D* as obtained 
from Eq. (34) together with contributions of light and charm quarks separately. We show 
the pion-mass dependence of the coupling constants and the charge radii in Fig. [3j 

In order to obtain the values of the observables at the chiral point, we perform linear and 
quadratic fits: 

/lin = aml + b, (35) 
/quad = aml + bml + c, (36) 

where a, b, c are the fit parameters. For all the observables we calculate here, both the 
linear and quadratic fits give consistent results, whereas the linear fit produces smaller 
errors. Our chiral-extrapolated result for gD*D-K, as can be seen in Table [1} is in good 
agreement with experiment. Table |TT] lists the chiral-extrapolated values of the meson cou- 
pling constants and charge radii. We observe that gD*D*p systematically lies above Qdop for 
all quark masses. According to SU(4) symmetry, we have goDp = gD*D*p which provides 
a good test for the amount of SU(4) breaking. We obtain {goDp/ gD*D*p)\m = 0.814(53) 
and {goDp/ gD*D*p)quad = 0.888(91) at the chiral point, which indicate an SU(4)-symmetry 
breaking around 20%. Comparing with gD*Dn, SU(4) breaking is much larger while SU(4) 
symmetry predicts goDp = gD*D*p = gD*Dn- Our results for the couphng constants can be 

12 




FIG. 3: dependence of the coupling constants goDp, gD*D*p and the charge radius (r^*) 

together with quark-sector contributions. The shaded regions show the chiral extrapolations with 
a linear form and their jackknife errors. 
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compared with those from QCD sum rules [281430] as goop = 2.9±0.4 and gD*D*p = 5.2±0.3, 
which imply a larger SU(4) breaking as compared to ours. Our value for gD*D*p is in agree- 
ment with that from QCD sum rules however there is a large discrepancy for goDp- 

It is expected from the quark model that the hyperfine interaction between the quark and 
the antiquark is repulsive in the vector mesons and attractive for the pseudoscalar mesons. 
We observe from Table [TT] that the charge radius of the D* meson is larger than that of D 
meson, indeed for all quark masses, in consistency with expectation from quark model. As 
it is seen in Fig. |3| the coupling constants and the charge radii of D and D* mesons tend 
to be coincident as the quark mass increases. This is a natural result because the hyperfine 
interaction, which is of the form -!I2-!l5_ jg reduced for a larger light-quark mass. We can 

' mqmq ' o o ^ 

also argue from this result that SU(4) breaking increases as the light-quark mass decreases. 
Indeed, similar relations hold for the p coupling constants in the heavy quark limit and one 
has goDp = gD*D*p (see Eq.(4.3) of Ref. [31]). 

Comparing the quark-sector contributions to the charge radii of D and D*, we find that 
the dominant contribution comes from the light quark. This implies that the large mass of 
the c quark suppresses the charge radii of D and D* mesons to smaller values compared to, 
e.g., charge radius of pion as (r^) = 0.452 fm^ [27]. This is also in qualitative agreement 
with the conclusion of Ref. [32j that the meson size decreases as the quark mass increases 
and that heavier quarks have distributions of smaller radius. The available literature on 
the electromagnetic properties of the D and D* mesons is limited [331 El]- Our result for 
the charge radius of the D meson is slightly below that from light-front quark model [31] as 
(rl) = 0.165;°:°i? fm^. 

A few comments on the systematic errors are in order now. To check the validity of 
the use of Clover action for the charm-quark in the case of coupling constants, we have 
repeated our measurements for k^. = 0.1216 and = 0.1232, which correspond to a change 
of ~ ±100 MeV in the charmonium mass. We have found that this leads to a change of 
less than 2% in the coupling constants as well as in the charge radii. Then our results are 
practically insensitive to a mild change in the charm-quark mass justifying the validity of 
Clover action in this case. 

For the finite- volume effect, the present spatial lattice extent is 32 units and the pion 
mass ranges from 0.136 to 0.322 in lattice unit, which gives 4.3 < rrij^L < 10.3. A rule of 
thumb is that serious effects seem to appear only when m^L < 4; therefore we expect small 
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TABLE II: The coupling constants goDp, 9D*D*p, the charge radius of D and D* mesons together 
with individual quark-sector contributions and meson masses {amp values are taken from [21j). We 
also give the D, D* and J/ip masses at different valence light quark masses. The chiral-extrapolated 
results are from linear and quadratic fits. 



K^al 9DDp {rh) (fm^) {rjj^^) (fm^) {rj^^^) (fm^) 


a mo aniji^ 


0.13700 5.57(27) 0.098(8) 0.226(20) 0.040(9) 
0.13727 5.03(28) 0.094(10) 0.232(29) 0.033(9) 
0.13754 5.30(30) 0.124(13) 0.350(50) 0.033(12) 
0.13770 4.85(41) 0.133(19) 0.308(85) 0.059(13) 


0.944(5) 1.453(5) 
0.919(4) 1.447(3) 
0.901(6) 1.434(6) 
0.896(10) 1.425(5) 


Lin. Fit 4.84(34) 0.138(13) 0.342(67) 0.051(11) 

Vc^uau. r It 4.yU(^0Dj U.iOZ(^/Dj U.o/Ul^lioj U.U(4(^iDj 




<'al 9D*D*p (rl,*)(fm2) (r2,. ,)(fm2) (r|,,^,) (fm^) 


a mo* amp 


0.13700 5.93(41) 0.106(12) 0.270(31) 0.035(13) 
0.13727 5.75(36) 0.113(17) 0.296(44) 0.036(14) 
0.13754 6.69(47) 0.167(25) 0.497(92) 0.044(21) 
0.13770 5.57(66) 0.169(30) 0.404(127) 0.075(26) 


1.006(7) 0.5060(30) 
0.981(8) 0.4566(36) 
0.971(7) 0.4108(31) 
0.940(9) 0.3895(94) 


Lin. Fit 5.94(56) 0.185(24) 0.475(94) 0.071(16) 
Quad. Fit 5.42(94) 0.192(43) 0.406(156) 0.096(29) 





finite-volume effects in our present calculations. 

We calculate only connected diagrams in this work. This means we neglect the effect 
of the disconnected diagram, where the external field is inserted in a sea-quark loop which 
in turn interacts with the meson two-point diagram with a gluon. The contribution of the 
disconnected diagrams is expected to cancel as long as we have an isovector field coupled 
to the meson as in the case of rho or pion couplings. A direct calculation of such diagrams 
appearing here only in the case of electromagnetic form factor, which has an isoscalar part 
as well as an isovector part, is difficult and numerically demanding. Their contribution has 
been found to be consistent with zero in the case of nucleon electric form factors [SSj- We 
expect them to give negligible contributions also in the case charm-meson electromagnetic 
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form factors. We also note that at higher-momentum transfers the effect of disconnected 
diagrams is further suppressed as a result of weaker coupling constant. 

IV. CONCLUSION 

Using the axial-vector coupling and the electromagnetic form factor of the D and D* 
mesons, we have computed the D*D7r, DDp and D*D*p coupling constants, which play 
an important role in describing the charm-hadron interactions in terms of meson-exchange 
models. We have also extracted the charge radii of D and D* mesons and determined 
the contributions of light and charm quarks separately. Our final results for the coupling 
constants as linearly extrapolated to the chiral point are 

gD'D. = 16.23(1.71), goDp = 4.84(34), gn*D'p = 5.94(56). (37) 

These suggest an SU(4) breaking around 20% for the p coupling constants, while it is more 
severely broken when they are compared to go'Dn- We have found that the SU(4) breaking 
gets larger as we decrease the light-quark mass. We have also calculated the charge radii of 
D and D* mesons and found that the dominant contribution to the charge radius comes from 
the light quark. We have found that the large mass of the c quark suppresses the charge 
radii of D and D* mesons to smaller values as compared to the charge radius of pion. 
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